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A  summary  and  the  details  of  the  research  performed  under  this 
grant  can  be  found  in  Attachment  1  (A  Summary  of  Recent  Topics  of 
Research  Results  on  Signal  Processing)  and  in  Attachment  2  (a 
listing  of  the  research  items) . 

Some  additional  work  is  detailed  in  the  following  reports  to  Rome 
Laboratory  which  were  all  sent  to  Dr.  Jon  Sjogren  who  was  the 
Program  Manager  for  this  grant. 

During  1992-93  the  principal  investigator  was  a  resident  at  the 
Rome  Laboratory  (URRP  Program)  and  the  result  of  the  work  is  the 
Reports  to  the  Rome  Laboratory  listed  in  Attachment  2. 


Attachment  1 


A  Report  to  Rome  Laboratory 
(1992-4a) 

A  SUMMARY  OF  RECENT  TOPICS 
OF  RESEARCH  RESULTS  ON 
SIGNAL  PROCESSING 


by 

E.  Emre 


Supported  by  the  AFOSR 


A  Report  to  Rome  Laboratory 
(1992-4a) 

A  SUMMARY  OF  RECENT  TOPICS 
OF  RESEARCH  RESULTS  ON 
SIOIAL  PROCESSING 


by 

E.  Smre 


Supported  by  AFOSR 


0.  SUMMARY  OF  TOPICS  FOR  ARRAY 
PROCESSING 

1.  TARGET  DENSITY  FUNCTION  ESTIMATION  VIA 

DIRECT  TECHNIQUES  (VIA  STANDARD  SIGNALS  AND 

ARRAYS) 

Provides  an  image  of  the  target  area  in 

terms  of  the  reflectivity  of  each  point  to  a 

physically  monochromatic  (or  very  narrowband) 
signal,  such  as  a  continuous  sine  wave,  or  a 
train  of  constant  amplitude  pulses  modulating 
a  sine  wave  (carrier) .  Hence  target 
identification/detection  is  possible  via  analysis 


of  this  image . 


TWO  lOklH  TYPES: 


la.  Time-invariant  (with  or  without  doppler 

shifts) .  Doppler  shift  models  the  velocity  of 

targets  assumed  to  move  with  constant  velocity. 

lb .  Time-varying : 

As  targets  move,  the  reflectivity  of  each 

point  in  space  varies  with  time. 

Estimation  of  a  time-varying  target 

density  function  provides 

i)  A  movie  picture  of  the  target  area  by 

as  many  snapshots  per  time  interval  as  desired 

ii)  Multitarget  Tracking  at  the  signal 
processing  level.  Replaces  Adaptive  Beamforming 

for  Tracking 
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2.  TARGET  DENSITY  FUNCTION  ESTIMATION 
AND  ANALYSIS  VIA  WAVELET 
AND  GABOR  THEORY 

Provides  multiresolution  spectral  analysis 
of  the  target  density  function  (image) . 

Also,  provides  new  techniques  for 

beamforming. 


3 .  TARGET  REFLECTIVITY  FREQUENCY 
RESPONSE  (TRANSFER  FUNCTION) 

AND/OR  IMPULSE  RESPONSE  ESTIMATION 
AT  EACH  ANGLE  SIMULTANEOUSLY. 

These  quantities  are  well  known  to  be 

target  signatures. 

At  a  fixed  range  (determined  via  such  as 
pulse  timing)  one  can  determine  these  quantities 
at  each  angle  simultaneously,  being  able  to 
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identify  each  target  sinultaneously. 


4 .  ANGLE-FREQUENCY  DENSITY 
FUNCTION  ESTIMATION 

These  techniques  enable  passive  sensors  to 
sort  the  incoming  signals  at  each  angle  and 
frequency  continually - 

Techniques  are  also  developed  to  determine 

dominant  (point)  angles,  frequencies  as 

extensions  of  MUSIC-PISARENKO  type  techniques. 

1.  TARGET  DENSITY  FUNCTION  ESTIMATION  VIA 
DIRECT  TECHNIQUES 

Consider  the  space  in  terms  of  (range=R, 
^“^ii^fictional  cosine)  coordinates. 
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Suppose  the  transmitted  signal  is 
j«ct 

W(t)  =  e 

or 

j«ct 

W(t)  =  e  p(t)  , 

where 

p(t)  *  a  train  of  pulses  of  constant 

amplitude. 

The  reflection  to  the  phase  center  from  a 

very  small  area  characterized  by  (R,R+ar)  and 

is  approximately 

n(R,p)W(t-aE^  (Ra/3)  (AR) 
c 

(c  *  speed  of  the  waves) 

The  set  of  reflectivity  constants 

for  each  R,/3  provides  an  image  of  the  target 

area.  As  ARr»0,  a/9-)0  the  total  reflected  signal 
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to  the  phase  center 


IS 


2R 


1  Rl 

(1)  y(0,t)  =  /  j.  h(R,^)W(t---)  RdRd^ 

-10  c 


where 


h(R,^)  =  lim  h(R,^) , 
AIHO 
A/3-iO 


and 


Rl  -  maximum  range  from  which  there  is  a 

reflection. 


For  simplicity  any  known  function  of  R 

such  as  due  to  medium,  integral  multiplier, 
etc.  will  be  included  as  a  factor  of 

h(R,^),  and  we  will  revvite 

(1)  as 

1  2r 

(2)  y(o,t)  =  /  /  h(R,iS)W(t- dRd)9. 

Target  Density  Function 
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The  sensor  at  point  x  receives 


1  2R  iSx 

(3)  y(x,t)  ^  S  I  h(R,jS)W(t - )  dRd^. 

-10  c  c 

Thus,  the  function 

h(R,^)  =  Target  Density  Function^  provides  an 
image  of  the  Target  area. 

If  the  target  area  changes  with  time, 

then,  the  reflectivity  of  each  point  (i8,R) 

varies  with  time.  This  time-varying  reflectivity 

is  represented  by  the 

h(t;/3,R)  =  Time-varvina  target  density  function. 

Thus,  if  h(t,/?,R)  can  be  estimated  for 

each  t,  or  for  as  many  values  of  t  as 

desired  per  processing  interval,  a  movie 

picture  of  the  target  area  is  obtained  with 
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«•  aany  Masar«Mnt8  for  each  target  point  per 


time  interval  as  desired.  Thus,  this  provides 

MULTI-TARGET  TRACKING  with  no  need  for  a 

posteriori  DATA-ASSOCIATION. 

Also,  provides  an  alternative  to  Adaptive 

Beamforming  for  tracking. 

An  alternative  (restrictive)  way  of 

modeling  time-variance  is  in  terms  of  doppler 

shift.  This  gives  rise  to  range-angle-doppler 

target  density  function 

h(u,R,/3)  with  «  =  a  particular  doppler  frequency 

corresponding  to  a  target  at 
(R,fi)  whose  velocity  is 
constant  and  yields  u  as  a 
doppler  shift. 
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COMPARISON  TO  PREVIOUS  MORE  ON 
TARGET  DENSITY  FUNCTION  ESTIMATION 

i)  Mostly  considers  density  functions  in 
range-doppler  coordinates 
h(u,R) ,  only. 

ii)  Most  such  techniques  require 

transmission  of  many  signals  (to  form 
a  basis) 

iii)  The  concept  of  a  truly  time-varying 
density  function  appears  in  range 
coordinates,  h(t,R) ;  with  no  general 
estimation  technique  given. 

We  have  developed  techniques  to  estimate 

target  density  functions  viewed  as  totally 

unknown  apriori 

i)  In  Range-Angle-Doppler  coordinates, 

h  (u ,  ^ ,  R) , 
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il)  Using  one  standard  radar  signal  and 
an  array . 

iii)  Also,  we  have  developed  techniques  to 
estimate  the  truly  time-varying  target 
density  functions 


h(t,/9,R)  . 

FURTHERMORE,  NOISE  IS  TAKEN  INTO  FULL 

CONSIDERATION  SEPARATELY. 


We  have  developed  techniques  for  estimation 


in  the  presence  of  noise  (e.g.,  receiver 


noise) .  Our  target  density  function  estimation 


techniques  utilize  Fourier  series  coefficients 


of  the  reflected  signal.  When  there  is  noise. 


we  obtain  the  MINIMUM  VARIANCE  estimate  of 


these  coefficients  (can  be  done  recursively  and 
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adaptively  using  a  Kalman  filter) . 

If  the  S/N  ratio  is  already  high,  one 


pulse  is  sufficient  for  target  density 


function  estimation. 
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2.  UKVSLBT  AMD  GABOR  THTORY  FOR 
KSTIMATKMI  AMD  ANALYSIS  OF  TARGET 
DENSITY  FONCTIOMS  AND  FOR  BEAMFORMING 


A  slight  modification  of  the  usual 

wavelet  theory  5  s  more  suited  to  array 
processing. 

A  wavelet  W(t)  is  a  function  where  its 

Fourier  transform  W(«)  satisfies 
oo  ) 

Y  j-  /  - ^ - d»*«o. 

*  —00  I  ***  i 

If  h(^)  is  the  angular  target  density 

function  at  a  fixed  range  R  (assumed  to  be  0 

for  simplicity),  and  if  W(t)  is  transmitted, 

the  sensor  at  point  x  receives 
1 

y(x,t)  =  /  g(/i)W(t-/3x)  d/3. 

-1 
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(take  c=l  for  simplicity) 

Since  g(/3)=o  for  ^^(-1,1), 

00 

y(x,t)  =  /  g(^)W(t-^x)  dfi 
-00 

=  (MODIFIED)  WAVELET  TRANSFORM  OF  g(/5)  . 

It  can  be  shovm  that 

00  00 

- /  /  y(x,t)W*(t-^x)dxdt 

"00 

=  INVERSE  TRANSFORM 

This  can  be  applied  to 

h(R,i8)  =  joint  range-angle  target  density  function 

with  respect  to  fi. 

Y(X,«)  =]^{y(X,t)}  =  W(o))H(««x) 

where 

H(«)  =  f{h(/9)}. 


Thus,  sensor  coordinate  x  is  the  dilation 


paraMter  for  nultirefcolution  spectral  analysis 


of  h(/9). 


SOME  OTHER  RESULTS 


i)  h(R,/9)  can  be  estimated  as  a  Fourier 


series  with  respect  to  a  wavelet 


ii)  very  easy  to  perform  beamforming. 


iii)  it  is  enough  to  have  finite  aperture 


(for  the  array),  and  finite  observation 


time.  Hence,  totally  practical. 


One  can  also  use  Gabor  theory  to 


estimate  h(/?,R)  and  spectral  analysis  with 


respect  to  R.  Many  results  (like  the  ed>ove 


for  wavelet  theory)  are  possible. 
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3 .  TARGET  REFLECTIVITY  IHPULSE  RESPONSE 
AND/OR  TRANSFER  FUNCTION  (SPECTRUM) 
ESTIMATION  AT  EACH  ANGLE  SIMULTANEOUSLY 

It  is  known  that  target  reflectivity  is 

frequency  dependent.  For  a  fixed  range,  if 

R(  «,/3)  is  the  frequency  dependent  reflectivity 

at  the  angle  /3,  techniques  are  developed  to 

estimate  R(  as  a  function  of 

H(  “V/3)  is  the  signature  of  the  target  at 

angle  fi.  Also,  technic[ues  are  developed  to 

estimate 


=  reflectivity  time  response  of  the 
target  at  angle  /3  to  an  impulse. 


Clearly,  in  linear  mode, 
f  (r(t,^)). 
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Itodiua  transfer  function  is  also  taken 

into  consideration. 

One  can  either  estimate  r(t,/9)  or  R(m,/3) 

directly,  for  each  /9  siumultaneously. 
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4.  DIRECTION-FREQUENCY  DECCNfPOSITION  OF 
SIGNAIf  RECEIVED  BY  AN  ARRY 
(PASSIVE  SENSORS) 

Two  main  types  of  results  are  obtained. 

i)  One  can  express  the  incoming  signal  to 
sensor  at  the  point  x  as 

1  00  j«  (t-^ 

y(x,t)  =  /  /  a(/9,u)e  ded/3. 

-1  -00 

.  Thus 

a(^,«)  is  the  direction-frequency  density 

function. 

Techniques  are  developed  to  estimate  aO,u} 

for  all 

ii)  One  may  be  interested  in  dominant 

directions  fi  for  a  given  « ,  as  well 
as  dominant  frequencies  u  for  a  given 
direction  p . 
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3?iK3hnlq[ues  are  developed  t:o  estimate 


dominant  directions  and  frequencies,  as 
extensions  of  MUSIC-PISARENKO  type  techniques. 

NOTE:  i)  is  shown  to  be  dual  to  the 

problem  of  wave-beamform  design:  i.e.,  find 

f(t,x)  =  the  signal  to  emit  from  the  sensor  at 

point  X,  such  that  in  each  direction  a 

possibly  different  signal  is  sent  from  the 
array,  i.e.,  the  points  at  the  angle  fi  at 

time  t  receive  the  wave-beam  form 

W(t,/9)  =  either  a  prespecified  desired 
function,  or  an  approximation 
to  it. 
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Attachment  2 


Report  to  Texas  Instrxments  1992  Series  a: 

2a)  A  New  Approach  to  Arrays  of  Sensors 

3a)  Simultaneous  Wave-Beam  form  Design  and  Angle-Frequency 
Density  Function  Estimation  (Temporal)  And  Spatial 
Decomposition 


Report  to  Texas  Instruments  1992  Series  b: 

lb)  Estimation  Of  Target  Density  Functions  In  Space  (Range- 
Angle)-  Doppler  Coordinates 

2b)  Range-Direction-Doppler  Target  Density  Fxmetion  Estimation 
By  Passive  Beamforming 

3b)  Joint  Range-Angle  Beamforming  With  Application  To  Estimation 
Of  Target  Density  Functions 

4b)  Estimation  Of  Time-Varying  Range-  Direction  Target  Density 
Functions  In  Nonstationary  Environments 


Report  to  Texas  Instiximents  Series  c: 

Ic)  Estimation  and  Multiresolution  Analysis  Of  Angular  and  Range- 
Angle  Target  Density  Functions/  Beamforming  Via  Wavelet 
Theory 

2c)  Estimation  And  Analysis  Of  Range  (Range-Angle)  Target  Density 
Functions  Via  (A  Modified)  Gabor  Theory 


Report  To  Rome  Laboratory  1992 


1)  Signal  Estimation  Techniques  For  Target  Density  Function 
Estimation  In  The  Presence  Of  Noise 

2)  A  Technique  For  Estimation  Of  Fourier  Series  Coefficients 
(And  Samples)  Of  A  Signal  From  It's  Noisy  Samples  Forming  The 
Usual  FFT  Bins 

3)  Virtual  Sensors  (Samples)  With  Application  To  Beamforming 
(Frequency  Filtering)  And  Approximation  Of  A  Continuously 
Distributed  Sensor  (Continuous  Function)  By  Point  Sensors 
(Time  Samples) 

4)  Remarks  On  Standard  (Physical)  Radar  Signals  And  Their  Fourier 
Series  Representation  For  Target  Density  Function  Estimation 

5)  A  Summary  Of  Work  On  Target  Density  Function  Estimation 

6)  Physical  Existence  Of  Target  Density  Functions 

7)  Existence  Of  Time-Varying  Target  Density  Functions  In 
Nonstationary  (Time-Varying)  Environments 

8)  Estimation  Of  Angular  Target  Reflectivity  Spectra  (Transfer 
Functions) 

9)  A  Direct  Technique  For  Estimation  Of  Angular  Target 
Reflectivity  Impulse  Response  Density  Function 

10)  Range-Angle  Target  Density  Function  Estimation  Via  A  Narrow 
Band  Signal 

11)  A  New  Technique  For  Estimation  Of  Time-Varying  Target  Density 
Functions  In  Nonstationary (Time-Varying)  Environments 


Reports  to  Rome  Laboratory  1992  Series  a: 

la)  Estimation  And  Analysis  Of  Target  Density  Functions  Via 
Wavelet  And  Gabor  Theory 

2a)  A  Svunmaxry  Of  Two  Techniques  For  Time-Varying  Target  Density 
Function  Estimation 

3a)  Advantages  Of  (Time-Varying)  Target  Density  Functions  For 
Radar  Signal  Processing 

4a)  A  Summary  Of  Recent  Topics  Of  Research  Results  On  Signal 
Processing 


Rome  Reports  1993 


1)  A  Modification  To  Certain  Target  Density  Function  Estimation 
Techniques 

la)  A  Svimmary  Of  Certain  Extrapolation  Techniques  In  The  Presence 
Of  Noise 


